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Abstract 

The method of complex characteristics of 
Garabedlan and Korn has been successfully used to 
design Ehockless cascades with solidities of up to 
one. A new code has been developed using this 
method and a new hodograph transformation of the 
flow onto an ellipse. This new code allows the 
design of cascades with solidities of up to two 
and larger turning angles, The equations of 
potential flow are solved In a complex hodograph- 
like domain by setting a characteristic Initial 
value problem and Integrating along suitable 
paths. The topology that the new mapping Intro- 
duces permits a simpler construction of these 
paths of Integration. 


The heuristic Idea behind the elliptic mapping 
Is the following. When the flow Is mapped Onto a 
circle, the upstream and downstream points of 
Infinity are mapped Into two Interior points of 
the circle, where two logarithmic singularities of 
the solution, a source and a sink, are located. 

By Increasing the separation between these two 
points the gap-to-chord ratio can be reduced. The 
limiting case of Infinite gap-to-chord ratio, or 
the Isolated airfoil, corresponds to the case In 
which both singularities coalesce. By mapping the 
flow onto an ellipse, the two singularities can be 
more widely separated, further reducing the gap- 
to-chord ratio. A parameter related to the 
eccentricity of the ellipse will be Introduced to 
control the solidity of the cascade. 


1. Introduction 


The design of supercritical blades for turbo- 
machinery relies heavily on the capability of 
producing a wide variety of two-dimensional blade 
sections, spanning from hub to tip and varying 
from upstream to downstream stages. A code to 
design these sections will have to cover a broad 
spentrum of design flow conditions. It should 
be able to handle both low and high solidities, 
chord-to-gap ratios gf up to two., as well as high 
Incidence angles and Mach numbers. 

The method of complex characteristics and 
hodograph transformation developed by Bauer, 
Garabedlan and Korn, [1], to design supercritical 
wing sections and cascades has been widely used. 
With this method, excellent shockless airfoils 
have been designed and tested In this country and 
abroad, In its latest version, the method solves 
the problem of finding a shockless airfoil with a 
given pressure distribution. 


With this new method we have achieved high 
solidities. A good resolution of points defining 
the body Is obtained both at the leading and 
trailing edge regions and larger Incldc.nce angles, 
with thinner airfoils, can be achieved. In the 
next section we review the method of, complex 
characteristics, as applied to our problem. 

Professor Paul R. Garabedlan of the Courant 
Institute of Mathematical Sciences proposed this 
problem and madn many suggestions during the 
development of the code, 

2. Hodograph Complex Equations 


The equations of potential flow In hodograph 
coordinates can be described by the Chaplygin 
system 



( 1 ) 


Because of the conformal mapping used by 
Bauer, Garabedlan and Korn to transform the hodo- 
graph domain of the flow onto a circle, a restric- 
tion exists In their method which does not allow 
the design of cascades with high solidities. 

Also, a poor resolution at the leading and trail- 
ing edge regions may arise, In cases of high 
Incidence angles, as a consequence of this 
mapping, 

In this paper a new design technique Is 
described, based on the same complex characteris- 
tics method, which uses a new mapping of the hodo- 
graph domain onto an ellipse. This new mapping 
leads naturally to the use of Tchebicheff poly- 
nomials, rather than trigonometric polynomials, 
to construct the solution of the equations. The 
equations of flow are Integrated along a new set 
of paths which are more In accord with the topol- 
ogy that the elliptic mapping introduces. 


♦Visiting Research Scientist presently at NASA 
Lewis Research Center, Fluid Mechanics and Acous- 
tics Division, Computational Fluid Mechanics 
Branch, Cleveland, Ohio 44135. 




In which f and ♦ are the potential, and 
stream functions, M Is the local Mach number, q 
and e are respectively the modulus and argu- 
ment of the velocity vector, and n Is the den- 
sity. This system of equations is of mixed type, 
elliptic and hyperbolic. In the transonic regime. 

In the method of complex characteristics, the 
variables are analytically extended Into the four 
dimensional domain of two complex variables, 
where the characteristic equations can always be 
solved. The system (1) can be written In the 
canonical form, [1], 



( 2 ) 


Itie two Independent variables K and n are 
arbitrary complex analytic functions of the char- 
acteristic coordinates 


S " lOp h - 19 


t - log h + 10 , 
where h 1s defined by 

^ VTT7/q' dq' 
h{q) » h(C*) e * 


^ , 

, ^ / Vl - M^/q* dq' (d) 

h(C*) m e 

and C* Is the rnltlcal speed. 

We recall that the coordinates s and t are 
conjuqated characteristics coordinates, [1], In 
the subsonic domain. This means that, for sub- 
sonic points, s and t are complex conjugate 
numbers when q and o i-e real. If the vari- 
ables c and n are defined In the form 


and the sonic locus. This part of the boundary 
will correspond to the subsonic part of the body. 
Points In the remainder of the domain 0 do not 
correspond to the real physical plane, and the 
real supersonic part of the body have to be found 
by searching for the zero stream line, 


The elliptic domain and the corresponding 
Tchebycheff polynomials can be defined In the 
following form. Consider the conformal mapping 

w - I 1 < iti < R (6) 

of a circular ring In the plane onto the ellipse 
with axes defined by the points 


and the slit (-1,1). We consider the class of 
analyfi functions In the ring which have a 
Laurent expansion of the form 



s - f(c). t - f(n), (5) 

with f analytic, then they are also conjugate 
coordinates. Under this transformation the sonic 
surface is mapped onto a cu'*ve of each char.c- 
terfstlc plane, called ihe ' anlc locus, [1]. 

In the hodograph i.jetaod, we are Interested In 
solving an Inverse orohiem. We construct a solu- 
tion with the corretL logarithmic singularities, 
which represent a sink and source at infinity, and 
we look at the zero stream line, on the phlsical 
real plane, as the possible candidate for the body 
which generates the flow. 


In which ag » a_g, for all positive n. The 
functions of the variable w 

... ( 8 ) 

are actually the Tchebycheff polynomials of the 
first kind, as can be easily verified- 

We observe that a function of the type (7) has 
the property 

F(e^») - F(e-''®), 0 < 9 < w ( 9 ) 


In the context of the complex system (2), the 
correct mithematical problem to solve Is the char- 
acterlstl<' Initial value problem where Initi.al 
data Is given 1n the two characteristic planes 
emanating from the Initial point (to.no) in 
the four-d1mensional space. The manipulation of 
this Initial data can and has lead In the past, 
[i), to the construction of shockless solutions. 

A new and fundamental approach In the application 
of hodograph methods was taken In [2], where the 
unknown domain of the flow is mapped onto the unit 
circle, by a conformal mapping of the type (5). 
This allows for a systematic prescription of the 
Initial characteristic data, based on a given 
pressure distribution over the airfoil. In the 
next section we describe how this Idea 1s 
implemented. 

3. Elliptic Conformal Transformation 

Following the Idea described In the last para- 
graph, we use the transformation (5) to map the 
flow onto an elliptic domain 0 of each coordi- 
nate characteristic plane. Because of this trans- 
formation, the coordinate characteristics are such 
that real subsonic po1nts_correspond to the domain 
of points of the form (e,e), when t Is a 
point enclosed by the portion of the boundary of 
D where t and n are conjugate coordinates 


This means that within this class of functions, we 
are Identifying points with the same real part In 
the positive and negative parts of the unit cir- 
cumference. Any analytic function on the ellipse 
In the (o-p1ane can then be written In the form 

■ £ ■ Z1 "n-" 

0 0 


This shows that the ellipse Is equivalent to the 
circular ring with the given identification of 
points on the unit circumference. This identifi- 
cation actually, makes unnecessary the slit Intro- 
duced by the mapping (6), making the elliptic 
liomeln singly connected. 

We adopt as our computational domain, 0, the 
circular ring 1 < Ki < R, with the previous 
Identification o7 poInTs. The parameter R will 
control the eccentricity of the ellipse and. In 
this way, the solidity of the cascade. Several 
advantages arise from this formulation, The most 
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Important one Is that we are now able to use Fast 
Fourier Transform to compute the coefficients of 
the mapping function 

fit) - 2 "n ('" * p:) 

of the flow ontv our computational domain 0. 

Other advantages of this formulation will be dis- 
cussed later. 

The supercritical design problem can be formu- 
lated as follows. Find ai* analytic (iolutlon 

»(C.n) ■ Rg{v^(t.n) log (n - n^) 


'!'(5,n) " Rg{*i{C.n) log (n - n^) 

♦git.n) log (n - ng) + ♦glt.n)} 


(IZ) 

of the system (2), with vi.fg.ea and ♦j.Vp.Vs 
regular functions, 1n the domain 0, where the 
mapping function f defined by (11), satisfies 
the Dirichlet condition 

Rg(f(5)j - log (h*(g)), ici •> R, (13) 


With the definitions given, h* and h become 
Identical In the subsonic part of Ipe domain D, 
where we take the parameter k equal to unity. 

In the remainder of this domain, h Is not real 
and the Introduction of the function h* becomes 
necessary, The consequence of using this 
Dirichlet condition Is that the designed airfoil 
will not achieve the proscribed pressure distribu- 
tion on the supersonic part of the airfoil. 

The boundary conditions (13) and (14) are non- 
llnearly coupled, The nonlinearity of the problem 
arises from the fact that the speed distribution 
Is given as a function of the arc length, and In 
(13) It Is needed as a function of the coordinate 
K on the boundary if >« R, The relationship 
between K and the arc length s, can only be 
determined when the potential function Is known. 

On the other hand, the solution of the boundary 
value problem (14) Is linked to the knowledge of 
the mapping function f. An Iterative procedure 
Is then needed to solve the problem. 

This Iterative procedure Is established by 
first computing an Incompressible solution 
f(e).'l'(t) on the elliptic domain. In terms 
of classic Jacobi elliptic functions. The In- 
compressible solution is conformal Invariant. 

We can then, establish the relation between the 
arc length s and the boundary variable t » Re'*, 
as we knov »(s). This allows computation of 
the roappiny function f, using F.F.T. to solve the 
Dirichlet problem (13) with prescribed values of 
h* at a number N of equidistant nodes on the 
circumference of radius R. 


and the stream function satisfies the boundary 
condition 


Rg{V(t,?)) » 0. Ki « R. 


(14) 


The stream function Is real valued for subsonic 
points and the boundary condition (14) is wr'1 
defined when a proper branch of the solution is 
taKen, The function h* is defined by 


h*(q) 



dq'/q' 


k > 0. 


Once the mapping function Is known, the system 
(2) can be Integrated, as t* are, then, known 
functions of the variables t and n. We re- 
mark that at each Iteration an analytic solution 
IS found. If the zero stream line Is not self 
Intersecting and has the proper closure at the 
trailing edge, It represents a shockless airfoil. 

4. Numerical Solution 


We concentrate In this section In the problem 
of finding a solution of the type (12) to the sys- 
tem (2), once the mapping function f has been 
determined. We introduce the notation 


(15) 

The regular part of the solution, can be 
expanded as a linear combination of particular 
solutions obtained by tak\ng as initial character- 
Isltic data a complete set of ortonormal func- 
tions, namely the Tchebychefr polynomials. The 
Coefficients of this linear coi«b1nation are then 
determined by imposing the boundary condition 
(14). This boundary value problem has been proved 
to be well posed when the system (1) Is reduced to 
Tricomi equation, [4], that Is, in the case of the 
transonic small disturbance equation. In the 
general case of the full potential equation that 
we are dealing with here, a low condition number 
for the matrix associated with the linear system 
In question, both when the domain 0 is a circle 
or an ellipse. Indicates the well posedness of the 
problem. 


L„(f,*) = v„ - 

The requirement that the functions (12) are 
solutions of the system (2) with the pairs 
(^^,i(i ^)1 B 1,3 being regular functioi^is leads 
to the homogeneous system, for each value 1 * 1,2 
of the Index 1, 


with 

on n = 1 “ 1.2 (17) 
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and t<) the Inhomogeneous system 



(ng*?) ■ Tjj(«) , n « 1 , 


TfiT 


ICO) 


(10) and where the Tn(u) are the Tchebicheff poly- 
nomlals described 1n (8), Once each of these 
Initial characterisfilc problems is nuraerlc'/rly 
solved, the coefficients b„ on (21) can be 
determined by Imposinq the boundary condition (14) 


The problem (16) with the condition (17) Is 
known as the Riemann probijm for that system. The 
solutions are the Riemann functions related to 
that pioblera and can be expllcity Integrated along 
the characterstic n k ni, In the form 


♦ |(c.ni) 


|/T+(c,n,) 


(19) 


and 1 If 1 , 2 . 

The complex constants Cj and C 2 remain 
at our disposal. These two sets of Riemann func- 
tions can be Integrated with a finite difference 
scheme, that will be described later. We can 
Impose on the t » C 1 , 1 ■ 1.2, characteristics 
the Initial data 


♦ ^(c^.n) * 1^^('n,t^) 


( 20 ) 


»^(c^,n) » v^(n.e,j). 


Once these functions are determined, we can 
Integrate the Inhomogeneous system (18). For 
that, we expand the solution In the form 


»3 




1 




bn* 


n’’3 


(n) 


( 21 ) 


where (» 3 °^ .♦ 3 °^ Is a solution of the Inhomogeneous 

problem (18) with the homogenous Initial character- 
tlstlc data 


♦5‘’)(t.n3) - (03,^ = 0, 


( 22 ) 


assigned on ^he two characteristics emanating from 
a point (C 3 .T 3 ), wh1ch_1s subsonic and far from 
the singularities (Cj,Ci) and Each pair 

of functions (f|^) >♦ 3 ^^) solutions of the 

homogeneous system (16) , with the Initial charac- 
teristic data 


A rectangular grid Is used to numerically 
solve each Initial characteristic problem. Two 
sides of this rectangular grid are formed by two 
paths of Integration, one In each characteristic 
plane, These paths allow one to reach the desired 
part of the flow to be computed. Three kinds of 
paths of Integration are used, subsonic, transonic 
and supersonic paths. 


To reduce the amount of computation, we divide 
the e.xterlor circumference Into eight sectorr. and 
we alviays take as initial point (Cg.cs) the 
Identical points (0,1) and ( 0 ,- 1 ) . From each 
singularity (ti.cj) and (C 2 ,t 2 )_we lay a path which 
goes directly to the polnt^c^.cg) to compute the 
Riemann solutions. The corresponding n~path 
will be the complex conjugate of the t-path. 

Once we rejich the characteristics through the 
point ( 53 , ? 3 ) we continue computing the 
Riemann solutions plus the solution of the in- 
homogeneous problem (18) In the way described. 


In those sectors In which all the nudes are 
subsonic points, we use subsonic paths. In this 
case the n-path is the complex conjugate of 
t-path. As t and n are conjugate coordi- 
nates at subsonic points, the diagonal of the 
rectangle corresponds to real subsonic points. 

The Initial data has previously been defined by 
the reflection laws (20), (22), and (23). In this 
form the computation can be reduced to the tri- 
angle below the diagonal, A path will consist, 
then, of 3 segment that starts at { 3 , ends at 
the outer circle, and continues with the circular 
arc corresponding to the sector In question. 


In those sectors, below the sonic locus, where 
all, or some, of the nodes do not correspond to 
real subsonic points we use transonic paths. The 
Idea Is that to reach those points beyond the 
sonic line, the two-dimensional manifold formed by 
the two paths has to avoid the sonic surface 
M(c,n) - 1 , where the equations become 111 
conditioned. So, In the case of those paths, we 
want the c-path and the n-path to not be con- 
jugates of each other. The Idea Is, as In [1], to 
traverse the corresponding sector In opposite 
directions for each path. 


Finally, the supersonic paths compute the real 
supersonic zone of the body. They were Introduced 
first by Swenson and they use the property that a 
point In the real supersonic zone can be reached 
by two characteristics starting at the sonic locus 
of each characteristic plane, [1], [5]. The com- 
putation for the supersonic zone Is done once the 
coefficients bn in ( 21 ) have been obtained. 

So the problem (18) is solved only once for these 
points. This makes affordable, both In terms of 
computing storage and C.P.U. time, the use of a 
much finer grid for the supersonic computation. 

The analytic solution is path independent, by the 
Cauchy Integral theorem, provided that we stay on 
the proper branch of the solution. 
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The finite difference scheme 

(24) 

’i,j ■ ’i,j-i " " ♦u-i) 

Is used to Integrate the equations. The average 
values t* are calculated with a predictor- 
corrector scheme which gives second-order accu- 
racy. A Richardson extrapolation to the zero 
limit Included In the code gives third-order 
accuracy. 

Two complex constants were left at pur dis- 
posal in the computation of the Rlemann func- 
tions. To determine them we require, first, the 
stream function to he single valued along the air- 
foil. This leads to the condition for the jump 
In ♦ 

C*) - . 0 125) 

On the other hand, the circulation r over the 
airfoil Is given. This Imposes the jump condition 
on f 

The location of the stagnation point and trailing 
edge Impose the two remaining conditions needed to 
determine the four real constants. 

With the potential and stream functions com- 
pletely determined, the body is calculated using 
the formula 

X + 1y .y ^ (d» + d*^ (27) 

By looking at the residuum of this function, when 
a loop 1s described around one singularity, say 
the source c = Kz, we obtain the repeating 'vetor 

ifig 

[x + 1y] - ~Zk ^ 

[*mf’2(^2’^2)) ^ ’m {*2(^2’^2))j‘ (2®) 

The modulus of this vector Is the gap, or distance 
between adjacent blades. If we consider. Instead, 
a loop which contains both singularities, we ob- 
tain the jump in the function x + iy around the 
airfnil 

[x + 1y] - -2« 


( 29 ) 


The real and Imaginary parts of this vector mea- 
sures the opening of the airfoil at the trailing 
edge In the x and y-d1rect1ons, respectively. 
This formula can be used to check the accuracy of 
the computation, when the actual coordinates of 
the body are obtained, 

S. Results 

A computer code has been written to Implement 
the method described. The purpose of the code 
1s to extend the possibilities of the Bauer- 
Garabedlan-Korn design code. We have followed, 
when possible, their method. A different strategy 
Is used In the numerical solution of the equa- 
tions, because of the different domain Into which 
the flow Is mapped. 

The code uses as Input a given speed distribu- 
tion, Other parameters used include the follow- 
ing: the external radius R of the circular 
ring, which controls the solidity, the number of 
nodes In the exterior circle, equal to the number 
of Tchebycheff polynomials used In the expansion 
of the regular part of the solution, the number of 
Iterations, which controls the total number of 
cycles (each o'sa Includes the calculation of the 
mapping function and a complete solution of the 
equations). Additionally, two grid sizes can be 
set, one for the subsonic-transonic computation, 
and the other, a finer one, for the supersonic 
computation. A Richardson extrapolation can be 
Included In the last Iteration. The Inlet Mach 
number Is specified In the case of compressor 
design, and the exit Mach number is specified In 
the case of turbines. 

A typical run takes 12 minutes of IBM 370-3033 
CPI) time. We remark that, while a blade design 
can Imply many runs to achieve the desired design 
conditions, most of these runs can be executed 
with a coarse grid, small number of nodes, and 
only one Iteration, which can reduce the C.P.U. 
time of each run to less than one minute, until 
the approximate design conditions are obtained. 

Figure 1 represents a high solidity cascade 
designed with the new code. A solidity of 2.05 
was obtained, with a turning angle of 42 degrees. 
The supersonic zone reached, though, has a mod- 
erate peak Mach number of 1.07. We use 128 *odes 
which delivers 193 computed points on the boty. 

It can be seen that a good resolution of points Is 
obtained both at the leading and trailing edge. 
Figure 2 represents a typical hodograph plane with 
the Integration paths. Figure 3 shows a highly 
staggered rotor tip section with an Inlet flow 
angle of 55 degrees. 

6. Conclusions 

A new design method has been developed for 
the design of supercritical cascades based on an 
elliptic confotmul transformation of the hodograph 
plane and the use of complex characteristics. 

With this new method we have been able to handle 
high solidity cascades. 

Because we use, as In the Bauer-Garabedlan- 
Korn method, an input pressure distribution the 
code can easily be coupled with a boundary layer 
calculation. In this way, the pressure distribu- 
tion can be modified until the desired separation 
criteria 1s met. 
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7. Refgrancei 


The effects of the strfe*"! tube convergence and 
radius change In the design of turbomachinery 
blades are well known. The Incorporation of those 
effects In the present code would enhance Its 
capabilities. A possible way to achieve this 
would be by reducing the three-dimensional problem 
to a set of two-dimensional problems by means of a 
Garlekln-type decomposition. The presence of a 
shroud In turbomachinery rows makes this method 
more appealing than In the case of Isolated 
airfoils. 
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Fig 1. High solidity stotot' blade 
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Fig 2, Hodograph plane 
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